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Determining the spectrum of hadronic excitations from Monte Carlo simulations requires the 
use of interpolating operators that couple to multi-particle states. Recent algorithmic advances 
have made the inclusion of multi-hadron operators in spectroscopy calculations a practical reality. 
In this talk, a procedure for constructing a set of multi-hadron interpolators that project onto 
the states of interest is described. To aid in the interpretation of simulation data, operators are 
designed to transform irreducibly under the lattice symmetry group. The identification of a set of 
optimal single-hadron interpolators for states with non-zero momenta is an essential intermediate 
step in this analysis. 
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Introduction 

The elucidation of the spectrum of hadronic excitations predicted by QCD is a long-term goal 
of lattice field theory. While a number of low-lying stable-particle masses are readily accessible, 
most of the low-lying spectrum has, until recently, lain beyond reach. The evaluation of correlation 
functions containing disconnected quark-line diagrams, which arise in the isosinglet-meson sector, 
is extremely challenging. Multi-particle correlators, needed to study states above threshold, are, in 
general, even more problematic. 

The stochastic LapH method, described in Ref. [[I]], is an efficient algorithm for evaluating 
hadronic correlators involving same-timeslice quark propagation. This algorithm facilitates the 
accurate evaluation of hadronic correlation functions using far fewer quark-matrix inversions than 
earlier methods. Moreover, using this approach, hadronic two-point functions factorize into source 
and sink contributions, greatly simplifying the evaluation of multi-hadron correlators in particular. 

The next step towards understanding the experimentally-accessible spectrum involves the cal- 
culation of the lowest-lying stationary-state energies on multiple lattice volumes. In a given sym- 
metry channel, multiple energy levels can be determined from the solutions of a generalized eigen- 
value problem involving matrices of hadronic two-point functions However, the success of this 
approach depends on the identification of appropriate interpolating operators. Considerable effort 
has already been invested in the construction of single-particle interpolators [||, Q|. In this note, we 
report on progress in constructing the multi-hadron interpolators needed above threshold. 

Although we are primarily interested in the spectrum of hadrons at rest, two or more hadrons 
in flight can combine to form a state with zero net momentum. Hence, as an intermediate step 
in our calculations, it is also necessary to construct interpolating operators for states with non- 
zero momenta. The analysis is further complicated by the fact that different flavor sectors become 
intertwined once multi-particle states are taken into account. For example, an isovector two-meson 
state might have one isosinglet constituent meson. To make sense of lattice data, it is essential to use 
interpolating operators with well-defined lattice quantum numbers. In particular, we use operators 
that transform irreducibly under the spatial symmetry group of the lattice. Our simulations are 
performed with 2 + 1 flavors of anisotropic clover fermions, and we are thus able to construct light- 
hadron operators with definite isospin or G-parity. Due to the finite temporal extent of the lattice, 
it is also advantageous to use meson operators that are irreducible under time reversal. Here, we 
focus on the spatial symmetries only. Incorporating the flavor and time-reversal symmetries is 
straightforward, and will be described in full in a forthcoming journal publication. 

Group theory 

In simulations on a cubic spatial lattice, the energies of zero-momentum states are commonly 
classified according to the irreducible representations (irreps) of the octahedral group, denoted Oh 
in Schonflies notation. The zero-momentum irreducible representations are a subset of the irreps of 
the space group of the cubic lattice, which is the semi-direct product of 0\ r and the group of lattice 
translations ST. The translation group is Abelian. Hence, all of its irreducible representations are 
one-dimensional. They are characterized by a lattice momentum p. Moreover, 5? is invariant under 
the full space group. Therefore, the method of induced representations is applicable, and the full 
set of space-group irreps can be deduced from the irreducible representations of the little groups of 
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the allowed lattice momenta. The little group of the lattice momentum p is the subgroup of lattice 
rotations which leaves p invariant. For mesons at rest, the relevant little group is the octahedral 
group. In the continuum, the little group of non-zero momenta is Coo V . C nv being the symmetry 
group of an n-sided regular pyramid. On the lattice, the little group depends on the orientation of 
the momentum with respect to the lattice axes. Hence, the little group of momenta along a lattice 
axis, (0,0,1), for example, is C^ v ; for planar-diagonal momenta, in the direction (0,1,1), say, the 
little group is Ci v \ and for momenta along a cubic diagonal, the little group is Cj, v . 

Periodic spatial boundary conditions are enforced in our simulations, and irreps for all lattice 
states with momenta satisfying |p| <4n/L can be induced from representations of the little groups 
listed above. The relevant little group for higher-momentum states moving in the (0, 1,2), (1,1,2) 
and equivalent directions is C2, the two-fold cyclic group. For the remainder of this report, we focus 
on states with momenta satisfying |p| < 4%/L. While we expect multi-particle states involving 
pions with higher momenta to play a role in the low-lying spectrum in current lattice simulations, 
the symmetry properties of these states are generally trivial to deduce. 

C$ v has five single- valued irreducible representations, which we label A\, A2, B\, B2 and E. 
We follow the Mulliken naming convention, so A and B denote one-dimensional representations, 
and E is two-dimensional. C% v has single-valued irreps A\,A%, B\, B2, and the irreps of C^ v are A\, 
A2 and E. The double-valued little-group irreps are one- or two-dimensional. To be consistent with 
the usual labeling convention for zero-momentum double- valued irreps [^J, we use G to denote 
two-dimensional representations, and F to denote one-dimensional fermionic irreps 1 . 

In the continuum and infinite-volume limits, non-zero momentum irreps are labeled by the 
absolute value of the helicity. The continuum representations are two-dimensional except for two 
zero-helicity irreps. The relationship between the various lattice little-group irreps and their con- 
tinuum counterparts and the irreps of Oh are determined by subduction. Table |] contains the sub- 
duction of the irreps of Oh to C^ v , and C^ v . The table tells us, for example, that a pion with 
momentum along a lattice axis transforms according to the A2 little-group irrep, and a vector me- 
son, in the T\ u irrep at rest, can appear in the A\ (scalar) or E irrep, depending on whether the 
momentum is parallel or perpendicular to the direction of polarization. 

Single-hadron operator construction 

Each space-group irrep is characterized by a set of lattice momenta closed under the action 
of Oj u and an irreducible representation of the corresponding little group. In practice, an operator 
basis for a space-group irrep is found by choosing a reference momentum p, and constructing a 
set of basis operators for the little group of p. The complete space-group basis is then obtained by 
applying a set of rotations which generate p*, the star of p, to this initial operator set. 

To construct little-irrep bases, we first identify gauge-invariant elemental operators with defi- 
nite momenta and flavor quantum numbers. The elementals incorporate different gauge-covariant 
displacements, but their transformation properties under lattice rotations are straightforward to de- 
duce. For example, the general expression for a meson elemental annihilation operator is 

*#(*,p) =I^ px (^J n)t ) (V) (W) (V), (i) 

'Our labeling convention differs from the convention in Ref. M, for example 
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Table 1: Subduction of the irreps of Oh onto the non-zero momentum little groups &x v , C2 V , and C^ v . 



where a (b) and A(B) are spin and flavor indices, and \f denotes a smeared quark field. The action 
of the gauge-covariant displacement operators is given by 

£>P (j>(x,t) = U k (x,t)---U k (x+(n- l)k, t) (x + nk, t) , (2) 

where U is a (stout-)smeared link variable. 

Under the action of the element R of the little group of p, (j) a (f,p) goes to 



M 
3=1 



The subset of N elementals {0} therefore generates a representation of the lattice little group, 
with matrices W, which is in general reducible. Basis operators for the constituent irreps are 
formed from linear combinations of these elementals as follows. First, given explicit representation 
matrices for the irrep A, F( A \ we define the N x N projection matrix 

ReG° 

where A labels a row of the A irrep and d\ is the dimension of the irrep. Uqd is the order of 
the group Gp 3 , the double cover of the little group of p. If the irrep A appears in the elemental 
representation W, and the elements of row a of P KX are not all zero, the linear superposition of 
elementals 



4 



Finite-momentum and multi-particle operators for hadron spectroscopy 



Justin Foley 









cm o 






ss 



SD 



DDL 



TDU 



TDO 



Figure 1: Representation of the spatial paths used in the meson elemental operators. All five displacement 
types are used in the construction of zero-momentum basis operators. Only single-site (SS) and singly- 
displaced (SD) operators are used in the non-zero momentum channels. 

transforms according to the row X of the irrep A. The rank of P is equal to the multiplicity of A 
in W. Starting with the operator for row A, the other operators in the little-inep basis are given by 



When constructing space-group operators, neither the reference momentum, p, nor the rota- 
tions used to generate p* are specified uniquely. However, a different choice of reference momen- 
tum or star rotations simply amounts to a change of basis, and is of no physical consequence. 

Single-hadron operator selection 

In order to project onto different hadronic excitations, a number of spatial-displacement com- 
binations are used in the operator sets. Fig. |T] shows the set of meson-elemental displacements 
currently in use. Of these, only the single-site and singly-displaced paths have been used in the 
non-zero momentum single-hadron operators to date. However, we are primarily interested in the 
lightest finite-momentum states, so we expect the simplest of the displaced interpolators to suffice. 

The operator-construction algorithm described above is readily automated, and one can easily 
generate a large number (in some of the zero-momentum irreps, hundreds or even thousands) of 
irreducible operators. It is therefore essential to apply a pruning procedure to the candidate operator 
sets to identify manageable subsets of clean operators that couple strongly to the lowest-lying states 
in each symmetry channel. The pruning consists of low-statistics analyses of the spectrum in each 
of the symmetry channels. First, the noisiest of the interpolators are identified and discarded. 
Cross-correlations between the remaining operators are examined in order to identify a subset of 
interpolators which couples well to a number of different states. 

Effective masses obtained from pruned operators in a non-zero momentum isoscalar meson 
channel can be seen in Fig. ||[ These measurements were performed on approximately 100 config- 
urations on a (1.9fm) 3 spatial lattice. These results correspond to the minimum allowed on-axis 
momenta. 

Multi-hadron states 

To form multi-hadron operators, representation matrices for the space-group operator bases 
are needed. Matrices for non-zero momentum irreps are obtained as follows. Having chosen a 
reference momentum p, and explicit bases for the little-group irreps {A} of p, with matrices r( A ), 




(6) 
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Figure 2: Lowest-lying effective masses in an isoscalar meson irrep containing states with the minimum 
on-axis lattice momenta. Effective masses were computed using pruned single-meson operators on approx- 
imately 100 configurations on a (1.9fm) 3 spatial lattice. The pion mass on these configurations is approxi- 
mately 400 MeV. Only single-site and singly-displaced operators were used. This particular lattice irrep is 
relevant for the rj meson. 



say, one identifies a subset of lattice rotations to generate p*, which we label R*. The components 
of the induced representation matrices are 

.(A) 



(a,pi),(j3,p2) 



0; p 2 + Rp 



where pi and P2 are momenta in p*, and /?, is the rotation in R* that takes the reference momentum 
p to p,. Hence if p2 = Rpi, R^ RR\ is a rotation in the little group of p. 

The procedure for constructing irreducible multi-hadron operators by combining irreducible 
single-particle operators is similar to the method used to construct single-hadron operators de- 
scribed previously. The direct products of single-particle irreps form (generally reducible) repre- 
sentations of the space group. These representations can be written in block-diagonal form, each 
block corresponding to a different relative orientation of the constituent single-particle momenta. 
To make two-particle operators for the rest spectrum, we need only consider the block with back-to- 
back single-hadron momenta, which itself forms a representation of O/,. Here, we use the notation 
[Ai,A2;p*] to denote this representation, where Ai and A2 are irreps of the little group of p*, the 
star of the single-particle momenta. The dimension of this two-particle representation is given by 
the product of the dimensions of Ai and A2, times the dimension of p*. 

Explicitly, the matrices corresponding to the zero-momentum block of the direct-product rep- 
resentations are 



< lA2;p,) (*) 



r (Ai; P *) , R) 



r (A2;p*) (R \ 

. («2,-p),(ft,-p') V ' 



(7) 
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where, on the left-hand side, we use a, /3 to denote the composite row and column indices (pt,\ , 052, p) 
and (j3i , /$2 , p') ■ In analogy to Eq. ^, we can define a projection matrix 

^ = ^~ I rfl{R)Mf^ (*)*. (8) 

Hence, if A appears in the decomposition of [Ai,A2;p*], and given the single-particle irreducible 
bases {C?( Al )}, jf?^ 2 '}, the composite operator 

P%0£\t, P >)0£\t,- P >) (9) 

(summation over j8, i.e., jSi, pi and p', implied) has zero net momentum, and transforms according 
to the row A of the octahedral-group irrep A. The multiplicity of the irrep A in the decomposition 
of [Ai,A2;p*] is given by the rank of the projection matrix P AX , which is simply the trace of P A ^. 
Averaging over all rows in A gives the multiplicity in terms of the characters of A and the little- 
group irreps: 

»a = — I \x {Al) (R)x (A2) (R)x iA> (R)\ ■ do) 

Using this formula, one finds, for example, that the six-dimensional representation [A2,A2; (0,0, 1)*], 
relevant for two pions with on-axis momenta, decomposes into A\ g ®E g © T\ u , and the twelve- 
dimensional pion-nucleon representation \A.2,G\\ (0,0, 1)*] contains the G\ g , G\ u , H g , H u irreps of 

o h . 

Work in progress 

To date, we have pruned single-hadron operators for all light-baryon, isovector-meson, isoscalar- 
meson, and kaon zero-momentum and non-zero momentum irreps. Two-particle coefficients for all 
zero-momentum irreps have been computed, and measurements of the single-hadron constituent 
components have been performed on lattices with estimated (3fm) 3 and (4fm) 3 spatial volumes 
and a pion mass of approximately 240 MeV. Results for the stationary-state spectra on these lat- 
tices will be presented in the near future. 

This work was supported by the U.S. NSF under awards PHY-0510020, PHY-0653315, PHY- 
0704171, PHY-0969863, and PHY-0970137, and through TeraGrid/XSEDE resources provided by 
the Pittsburgh Supercomputer Center, the Texas Advanced Computing Center, and the National 
Institute for Computational Sciences under grant numbers TG-PHY100027 and TG-MCA075017. 
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